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The temperature and electric field dependence of electronic contribution to the thermal conductivity (TC) 
of a granular superconductor is considered within a 3D model of inductive Josephson junction arrays. In addi- 
tion to a low-temperature maximum of zero-field TC k(T, 0) (controlled by mutual inductance Lo and normal 
state resistivity R n ), the model predicts two major effects in applied electric field: (i) decrease of the linear 
TC, and (ii) giant enhancement of the nonlinear (i.e. VT - dependent) TC with A/t(T, E)/k(T, 0) reaching 
500% for parallel electric fields E ~ Tt (Et ~ 5o|VT is an "intrinsic" thermoelectric field). A possiblity of 
experimental observation of the predicted effects in granular superconductors is discussed. 



PACS: 74.25.Fy, 74.50.+r, 74.80.Bj 

1. Introduction. Inspired by new possibilities 
offered by the cutting-edge nanotechnologies, the ex- 
perimental and theoretical physics of increasingly so- 
phisticated mesoscopic quantum devices (heavily based 
on Josephson junctions and their arrays) is becoming 
one of the most exciting and rapidly growing areas of 
modern science [[l], || Ej]- ^ n addition to the tradi- 
tional fields of expertise (such as granular superconduc- 
tors Q), Josephson junction arrays (JJAs) are actively 
used for testing principally novel ideas (like, e.g., topo- 
logically protected quantum bits |]) in a bid to solve 
probably one of the most challenging problems in quan- 
tum computing. Though traditionally, the main empha- 
sis in studying JJAs has been on their behavior in ap- 
plied magnetic fields, since recently a special attention 
has been given to the so-called electric field effects (FEs) 
in JJs and granular superconductors 0, [|, ||, 0, ||, 
The unusually strong FEs observed in bulk high-T c 
superconducting (HTS) ceramics [Q (including a sub- 
stantial enhancement of the critical current, reaching 
AI c {E)/I c (0) = 100% for E = 10 7 V/m) have been at- 
tributed to a crucial modification of the original weak- 
links structure under the influence of very strong electric 
fields. This hypothesis has been corroborated by further 
investigations, both experimental (through observation 
of the correlation between the critical current behavior 
and type of weak links || ) and theoretical (by studying 
the FEs in SNS-type structures || and <i-wave granular 
superconductors ffl). Among other interesting field in- 
duced effects, one can mention the FE-based Josephson 
transistor m and Josephson analog of the magnetoelec- 
tric effect |9|] (electric field generation of Josephson mag- 
netic moment in zero magnetic field). At the same time, 
very little is known about influence of electric fields on 



thermal transport properties of granular superconduc- 
tors. In an attempt to shed some light on this intresting 
and important (for potential applications) problem, in 
this Letter we present a theoretical study of the electric 
field and temperature dependence of electronic contri- 
bution to thermal conductivity (TC) n of a granular 
superconductor (described by a 3D model of inductive 
JJAs). As we shall see below, in addition to a low- 
temperature maximum of zero-field TC n(T, 0) (con- 
trolled by the mutual inductance Lq and normal state 
resistivity R n ), the model predicts unusually strong (gi- 
ant) field-induced effects in the behavior of nonlinear 
(i.e. VT - dependent) TC. In particular, the absolute 
values of the TC enhancement Ak(T, E) / k(T, 0) are es- 
timated to reach up to 500% for relatively low (in com- 
parison with the fields needed to observe a critical cur- 
rent enhancement j| |j| ) applied electric fields E match- 
ing an intrinsic thermoelectric field Et = iSo|VT|. The 
estimates of the model parameters suggest quite an op- 
timistic possibility to observe the predicted effects in 
granular superconductors and JJAs. 

2. The model. To adequately describe a thermo- 
dynamic behavior of a real granular superconductor for 
all temperatures and under a simultaneous influence of 
arbitrary electric field E and thermal gradient VT, we 
consider one of the numerous versions of the 3D JJAs 
models based on the following Hamiltonian 



n{t) =n T (t) +n L (t) + n E (t), 



where 



N 



n T (t) = ^2j ij [i-cos(pi j (t)] 



(i) 



(2) 
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is the well-known tunneling Hamiltonian, 



2L,, 



(3) 



accounts for a mutual inductance between grains 
(and controls the normal state value of the thermal con- 
ductivity, see below) with $>ij(t) — (h/2e)(f>ij(t) being 
the total magnetic flux through an array, and finally 



n E (t) = -p(t)E 



(4) 



describes electric field induced polarization contribu- 
tion, where the polarization operator 



N 



p(t) = -2eJ2ni(t>i 



(5) 



i=l 



Here nt is the pair number operator, and r; is the coor- 
dinate of the center of the grain. 

As usual, the tunneling Hamiltonian Tixit) describes 
a short-range interaction between iV superconducting 
grains, arranged in a 3D lattice with coordinates = 
(%i,yi,Zi). The grains are separated by insulating 
boundaries producing temperature dependent Joseph- 
son coupling J i:j (T) = J i:j (O)F(T) with 



W^tanh 



A(T) 



2k B T 



(6) 



and Jy(0) = [A(0)/2}(R /R lj ) where A(T) is the tem- 
perature dependent gap parameter, Rq = h/Ae 2 is the 
quantum resistance, and R t j is the resistance between 
grains in their normal state assumed [ flO[ to vary ex- 
ponentially with the distance between neighboring 
grains, i.e. = R~ l exp(— rij/d) (where d is of the 

order of an average grain size) . 

As is well-known [HJ, a constant electric field 
E and a thermal gradient VT applied to a JJA cause 
a time evolution of the initial phase difference = 
d>4 — d>H as follows 



(E, VT)f 



(7) 



evolution with a mutual inductance Lij. For simplic- 
ity, in what follows we limit our consideration to the 



longitudinal component of the total thermal flux Q(i) 
which is defined (in a q-space representation) via the 
total energy conservation law as follows 



Q(t) 



lim 



rHq(*) 



where ?i q = dTtq/dt with 

Wq(t) = ~ J d 3 xe itir H(r,t) 



(8) 



(9) 



is properly defined normaliza- 
we made a usual substitution 



Here v = 87tg? 3 
tion volume, and 
wT,ij A (nj,t) ±fd 3 xA(r,t) valid in the long- 
wavelength approximation (q — > 0). 

In turn, the above-introduced heat flux Q(t) is re- 
lated to the appropriate components of the linear ther- 
mal conductivity (LTC) tensor follows (hereafter 
{a, 13} = x,y,z) 



k q/3 (T,E) 



1 

V 



d< Q a > 



VT=0 



where 



< 



Qa> = - [ dt< Q a (t) > 
T JO 



(10) 



(11) 



Here V is a sample's volume, r is a characteristic 
Josephson time for the network, and < ... > denotes 
the thermodynamic averaging over the initial phase dif- 



ferences (jf-j 



< A(4>%) >=^f o X\d^, J A(^)e'^ (12) 



with an effective Hamiltonian 



T dt 

T 



d A X 



H(r,t) 



(13) 



Here, j3 = l/k B T, and Z = /* fTy - d^e' 13110 is the 
partition function. The above-defined averaging proce- 
dure allows us to study the temperature evolution of the 
system. 

Taking into account that in J JAs p"l[ oc R^j , we 



Here Wjj = 2e(E — E-r)rij//i where Et = SqVT is an obtain L 
" intrinsic" thermoelectric field with Sq being a zero-field 
value of the Seebeck coefficient. 

3. Linear thermal conductivity (Fourier law). 
We start our consideration by discussing the tempera- 
ture behavior of the conventional (that is linear) thermal 
conductivity of a granular superconductor in arbitrary 
applied electric field E paying a special attention to its 



y — Lq exp(rjj / d) for the explicit r-dependence 
of the weak-link inductance in our model. Finally, 
in view of Eqs.(l)-(13), and making use of the usual 



" phase- number" commutation relation, 



= i5. 



13 > 



we find the following analytical expression for the tem- 
perature and electric field dependence of the electronic 
contribution to linear thermal conductivity of a granular 
superconductor 

K af3 (T,E) = K O [8 a pr,{T,e) + f3 L {T)u{T,e)f a0 {e)) (14) 
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where 



fafi(e) = - [8 a pA(e) - e a epB(e) 



with 



and 



A(e) 



5 + 3e 2 
(1 + e 2 ) 2 



• tan 



B(e) 



3e 4 + 8e 2 



e 2 (l + e 2 ) 3 



tan 



(15) 



(16) 



(17) 



Here, k q = Nd 2 S $ /VL , I3 L {T) = 2irI c (T)L Q /<S> 
with I C (T) = (2e/h)J(T) being the critical current (we 
neglect a possible field dependence of I c because, as 
we shall see below, the characteristic fields where ther- 
mal conductivity exhibits most interesting behavior are 
much lower than those needed to produce a tangible 

change of the critical current |Q); e = ^/e 2 + e 2 , + e 2 

with e a = E a /Eo, and Eq = h/(2edr) is a character- 
istic electric field. In turn, the above- introduced "or- 
der parameters" of the system, r](T, e) =< > and 
v(T,e) =< sine/)?- >, are defined as follows 



TV 

2 



E 



7T ^ (2n+ l) 2 

n=0 ^ ; 



We) 



and 



where 



Pe(T,c) 



i/(T,e) 



0J(T) 



sinh Pe 

I3 e Iq(Pe) : 



1 



■ tan 



(18) 



(19) 



(20) 



Here J(T) = J(0)F(T) with J(0) = (A /2){R /R n ) 
and i ? (T) given by Eq.(6); I n (x) stand for the appro- 
priate modified Bessel functions. 

3.1. Zero-field effects. Turning to the discussion 
of the obtained results, we start with a more sim- 
ple zero-field case. The relevant parameters affecting 
the behavior of the LTC in this particular case in- 
clude the mutual inductance L and the normal state 
resistance between grains R n . For the temperature 
dependence of the Josephson energy (see Eq.(6)), we 
used the well-known |L2| approximation for the BCS 
gap parameter, valid for all temperatures, A(T) = 

A(0) tanh ^y^J^~^j with 7 = 2.2. 

Despite a rather simplified nature of our model, it 
seems to quite reasonably describe the behavior of the 
LTC for all temperatures. Indeed, in the absence of an 
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Fig.jl]. Temperature dependence of the zero-field linear 
thermal conductivity kl (T, 0) / n n for different values of 
the dimensionless parameter /3l(0). 



applied electric field (E — 0), the LTC is isotropic (as 
expected), K a p(T,0) — S a /3KL(T, 0) where «l(T, 0) = 
Ko[v(T,0) + 2/3l(T)i>(T, 0)] vanishes at zero tempera- 
ture and reaches a normal state value k„ = kl(T c , 0) = 
(tt/2)kq at T = T c . Figure |l| shows the temperature 
dependence of the normalized LTC k l {T, 0) / n n for dif- 
ferent values of the dimensionless parameter /3l(0) = 
2irl c (0)Lq/^o . As it is clearly seen, with increasing of 
this parameter, the LTC evolves from a flat-like pattern 
(for a relatively small values of L ) to a low-temperature 
maximum (for higher values of /3l(0)). Notice that the 
peak temperature T p is practically insensitive to the 
variation of inductance parameter Lq while being at the 
same time strongly influenced by resistivity R n . Indeed, 
the presented here curves correspond to the resistance 
ratio r n = Ro/R n = 1 (a highly resistive state). It can 
be shown that a different choice of r n leads to quite a 
tangible shifting of the maximum. Namely, the smaller 
is the normal resistance between grains R n (or the bet- 
ter is the quality of the sample) the higher is the tem- 
perature at which the peak is developed. As a matter of 
fact, the peak temperature T p is related to the so-called 
phase-locking temperature Tj (which marks the estab- 
lishment of phase coherence between the adjacent grains 
in the array and always lies below a single grain super- 
conducting temperature T c ) which is usually defined via 
an average (per grain) Josephson coupling energy as Jll| 
J(Tj,r n ) — ksTj. In particular, for T ~ T c , it can be 
shown analytically that Tj(r n ) indeed increases with r n 
as Tj(r n )/T c ~ r„/(l + r n ). 

3.2. Electric field effects. Turning to the discus- 
sion of the LTC behavior in applied electric field, let us 
demonstrate first of all its anisotropic nature. For sim- 
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Electric field dependence of the nonlinear ther- 
mal conductivity (T, E) /k*£ (T, 0) for different val- 
ues of the applied thermal gradient er = So]VT|/ Eo 
(ex = 0.2; 0.4; 0.6; 0.8; 1.0, increasing from bottom to 
top). Inset: Electric field dependence of the lin- 
ear thermal conductivity Kl(T, E)/ Kl(T,Q) for parallel 
(E[|VT) and perpendicular (E_LVT) configurations. 



plicity (but without losing generality), we assume that 
E = (£,0,0) and VT = (V X T, V y T, 0) . Such a choice 
of the external fields allows us to consider both paral- 
lel k xx (T,E) and perpendicular K yy (T,E) components 
of the LTC corresponding to the two most interesting 
configurations, E||VT and EJ_VT, respectively. Inset 
in Figure || demonstrates the predicted electric field de- 
pendence of the normalized LTC Kl(T, E)/k l (T 7 0) for 
both configurations taken at T = 0.2T C (with r„ = 1 and 
— 1)- First of all, we note that both components 
of the LTC are decreasing with increasing of the field 
E/Eq. And secondly, the normal component K yy de- 
creases more slowly than the parallel one K xx , suggesting 
thus some kind of anistropy in the system. In view of the 
structure of Eq.(14), the same behavior is also expected 
for the temperature dependence of the field-induced 
LTC, that is Ak l (T, E)/k l (T, 0) < for all fields and 
temperatures. In terms of the absolute values, for T = 
0.2T C and E = E , we obtain [Ak l (T, E)/k l (T, 0)} xx = 
90% and [AK L (T,E)/K L {T,0)} yy = 60% for attenuation 
of LTC in applied electric field. 

4. Nonlinear thermal conductivity: Giant 
field-induced effects. Let us turn now to the most 
intriguing part of this paper and consider a nonlin- 
ear generalization of the Fourier law and very unusual 
behavior of the resulting nonlinear thermal conductiv- 
ity (NLTC) under the influence of an applied electric 
field. In what follows, by the NLTC we understand 



a VT-dependent thermal conductivity K^g(T, E) 
n a js(T, E; VT) which is defined as follows 



S ^(T,E) = -- 



d< Q a > 



(21) 



vivo 



with < Q a > given by Eq.(ll). 

Repeating the same procedure as before, we obtain 
finally for the relevant components of the NLTC tensor 



K^(T,E) = K O [S a0 r,(T,e eff ) 

+{3 L (T)V(T, 6eff)D a(3 (6eff)}, 



(22) 



where 



D a p{eeff) = f a p{teff) + e T g af3l {t ef f) (23) 



ith 



ff Q /3 7 (e) = g[(<5 Q /3£ 7 + Sa^ep + 5~ (l3 e a )B(e) (24) 
+3e Q e /3 e 7 C(e)] 



and 



C(e) = 



lie 2 - lie 4 



3e b 



e 4 (l + e 2 ) 4 



• tan 



(25) 



Here, e"jj — e a — e T where e a — E a /E and e T — 
E^/Eq with Ej, — 5oV Q T; other field-dependent pa- 
rameters (?y, v 1 B and f a p) are the same as before but 
with e — > e e ff. 

As expected, in the limit Et — ► (or when E 3> 
Et), from Eq.(22) we recover all the results obtained 
in the previous section for the LTC. Let us see now 
what happens when the "intrinsic" thermoelectric field 
Ex = SoVT becomes comparable with an applied elec- 
tric field E. Figure || (main frame) depicts the result- 
ing electric field dependence of the parallel component 
of the NLTC tensor k* x l (T,E) for different values of 
the dimensionless parameter ct = Et/Eq (the other 
parameters are the same as before). As it is clearly 
seen from this picture, in a sharp contrast with the 
field behavior of the previously considered linear TC, 
its nonlinear analog evolves with the field quite dif- 
ferently. Namely, NLTC strongly increases for small 
electric fields (E < E m ), reaches a pronounced maxi- 
mum at E = E m = |-Et; and eventually declines at 
higher fields (E > E m ). Furthermore, as it directly 
follows from the very structure of Eq.(22), a similar 
" reentrant-like" behavior of the nonlinear thermal con- 
ductivity will occur in its temperature dependence as 
well. Even more remarkable is the absolute value of 
the field-induced enhancement. According to Figure || 
(main frame), it is easy to estimate that near maximum 
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(with E — E m and E T = E ) and for T = 0.2T C , one 
gets A K » x l (T,E)/k% x l (T,0) ~ 500%. 

5. Discussion. To understand the above-obtained 
rather unusual results, let us take a closer look at the 
field-induced behavior of the Josephson voltage in our 
system (see Eq.(7)). Clearly, strong heat conduction re- 
quires establishment of a quasi-stationary (that is nearly 
zero-voltage) regime within the array. In other words, 
the maximum of the thermal conductivity in applied 
electric field should correlate with a minimum of the to- 
tal voltage in the system, V(E) = (A) < >= 
V (e - ct) where e = E/E and V = E d = h/2er is a 
characteristic voltage. For linear TC (which is valid only 
for small thermal gradients with ct = Et/Eq <C 1), the 
average voltage through an array Vl{E) ~ Vq{E/E ) 
has a minimum at zero applied field (where LTC indeed 
has its maximum value, see the inset of Figure ||) while 
for nonlinear TC (with ej 1 — 1) we have to consider the 
total voltage V(E) which becomes minimal at E = Et 
(in a good agreement with the predictions for NLTC 
maximum which appears at E = §.Et> see the main 
frame of Figure |^) . 

To complete our study, let us estimate an order of 
magnitude of the main model parameters. Starting with 
applied electric fields E needed to observe the above- 
predicted nonlinear field effects in granular supercon- 
ductors, we notice that according to Figure |[ the most 
interesting behavior of NLTC takes place for E ~ Eq. 
Taking d ~ 10^m and r ~ 10~ 9 s for typical values of 
the average grain size and characteristic Josephson tun- 
neling time (valid for conventional JJs H] and HTS 
ceramics (l(J), we get E = h/(2edr) ~ 2 x 10~ 2 V/m 
for the characteristic electric field (which is surprisingly 
lower than the typical fields needed to observe a criti- 
cal current enhancement in HTS ceramics ||, ||). On 
the other hand, the maximum of NLTC occurs when 
this field nearly perfectly matches an "intrinsic" ther- 
moelectric field Et = Sb|VT| induced by an applied 
thermal gradient, that is when E ~ Eq ~ Et- Using 
So — Q.huV/K for the zero-field value of the linear See- 
beck coefficient 0, Q, we obtain \VT\ E ^ E /S ~ 
4 x lQ 4 K/m for the characteristic value of an applied 
thermal gradient. Finally, taking as an example granu- 
lar aluminum films with phonon dominated heat trans- 
port |l§ (with K ph (T) ~ 2 x 10- J W/mK at T = Tj ~ 
0.2T C ), let us estimate the absolute value of the pre- 
dicted here zero- field electronic contribution n e {T) = 
Kl(T, 0) at T = 0.2T C . Recalling that within our model 
the scattering of normal electrons is due to the pres- 
ence of mutual inductance between the adjacent grains 
Lq, and assuming that Q Lq ~ nod ~ 4ir x 10~ 12 H 
and V ~ Nd 2 l with I ~ 0.5mm (I is a film's thickness), 



we obtain n e (T = 0.2T C ) ~ f3 L (0) x IQ^W/mK for a 
rough estimate of the electronic contribution to the dis- 
cussed here inductance-driven effect. Correspondingly, 
we get K e (Q.2T c )/K, ph {0.2T c ) ~ /3 L (0)/2 for the ratio, 
where /3 L (0) = 2irI c (0)L /% (for example, (3 L (0) ~ 4 
for J c (0) = 10~ 4 A). Thus, depending mainly on the 
value of the critical current I c (0) and mutual inductance 
between adjacent grains Lq, the thermal conductivity of 
specially prepared granular alumina films will be dom- 
inated by either phonon (for small /3l(0)) or electronic 
(for large (3l(0)) contribution. Undoubtedly, the above 
estimates suggest quite a realistic possibility to observe 
the predicted non-trivial behavior of the thermal con- 
ductivity in granular superconductors and artificially 
prepared Josephson junction arrays. We hope that the 
presented here results will motivate further theoretical 
and experimental studies of this interesting problem. 
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